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Abstract. We derive new constraints on the zeros of Airy functions by using the 
so-called quantum bouncer system to evaluate quantum-mechanical sum rules and 
perform perturbation theory calculations for the Stark effect. Using commutation 
£^ ' and completeness relations, we show how to systematically evaluate sums of the form 

S p (n) — X)/t=4n V(Cfc — C«) P j f° r natural p > 1, where — £„ is the n th zero of Ai((). 
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1. Introduction 

The study of special functions and orthogonal polynomials as solutions of some of the 
most important and physically relevant model problems in quantum theory has formed 
one of the foundations of quantum mechanics since its earliest days. Familiar examples 
include Hermite and Laguerre polynomials for the harmonic oscillator and Coulomb 
problems, spherical harmonics (Legendre polynomials) for angular momentum, Bessel 
functions (in two and three dimensions) for both scattering and bound-state problems, 
Mathieu functions for the quantum-mechanical version of the pendulum, and even 
trigonometric functions for the infinite square well problem. 

The Airy function p] has found many applications in classical physics, especially 
in optics and fluid mechanics, and as an important tool in the derivation of the WKB 
approximation in quantum mechanics. Airy functions are also central to the 'quantum 
bouncer' problem, a point mass subject to the potential 

V(x) = { Fx l OIX> ° . (1) 

u |oo forx<0 v ; 

This problem is amenable to a variety of approximate treatments (WKB approximation, 
variational methods, numerical approaches including the shooting method), all of which 
can be compared to the exact solutions, given in terms of the Airy function Ai((). This 
system has also received considerable recent interest as it is the simplest model for the 
"quantum states of neutrons in the Earth's gravitational field" observed by Nesvizhevsky 
et al. [2]. 

Gea-Banacloche [3] investigated the time- dependence of localized wave packets in 
such a potential, motivated by the possible application to ultra-cold atoms dropped onto 
an 'atomic mirror,' and demonstrated the existence of quantum wave packet revivals [I] 
which have been observed in a wide variety of quantum-mechanical systems. He also 
found for the first time numerical evidence for expressions for the normalization of the 
'quantum bouncer' energy eigenstates and dipole matrix elements, (n\x\k), all of which 
were quickly shown to be analytically correct by Vallee [5] using earlier published work 
from the mathematical [6] and scientific literature [7] on Airy functions. Goodmanson 
[8] extended these results to produce closed-form expressions for the matrix elements 
of arbitrary powers of position, (n\x q \k), using recursion relations. The lowest-order 
(dipole and quadrupole) off-diagonal elements, {n\x 1,2 \k), have very simple forms and 
were found to depend on inverse powers of the combination (( n — ( k ) where — ( n is the 
n th zero of Ai((). 

Such differences are intrinsically related to the corresponding differences in energy 
eigenvalues for the stationary states, E n — E&, since for the 'quantum bouncer,' the 
quantized energy eigenvalues are given directly by 

(h 2 F 2 \ 1 ^ 

E n = ( n £ where £ = — — . (2) 

y 2m J 

The association of such energy differences in conjunction with dipole (and higher) matrix 
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elements for a quantum-mechanical system is most familiar from energy-weighted sum 
rules, such as the Thomas-Reiche-Kuhn (TRK) formula [9] 

h 2 

]T(£ fc - E n )\{n\x\k)\ 2 = — (TRK sum rule) . (3) 

k lm 

Because of the simple nature of the expressions found in Refs. [3] - [8], such sum rules 

automatically provide constraints on sums of inverse powers of differences of the zeros 

of Airy functions, and a systematic study of such relationships is the topic which we will 

examine here. (We note that Sukumar [10], [TTJ has found constraints on combinations 

of inverse powers of energy eigenvalues using Green's function techniques for several 

familiar model systems [10], including the linear potential involving the Airy function 

[TT] studied here. Our approach is different, however, and we find new constraints on the 

C, n which have not, to our knowledge, appeared in the mathematical physics literature.) 

In Section [2] we briefly review the quantum-mechanical solution of the 'quantum 

bouncer' problem in terms of Airy functions, and use a number of well-known quantum 

sum rules to find new constraints on Airy function zeros. Specifically, we find closed 

form expressions for the quantity 

SM - £ 77-^ <«) 

for numerous natural (positive integer) values of p. It has also been noted that energy- 
difference weighted sum rules have a mathematical structure similar to second-order 
perturbation theory [12] and we will also use results for the second-order Stark shift 
(the energy change due to the addition of an external constant field) in the same way. 
In Section [3] we then show how to systematically derive closed form expressions for the 
S p (n) in Eqn. (jlj) using commutation relations and closure methods (inserting a complete 
set of states) for all natural values of p > 1, and exhibit specific results for p — 2, ..., 11 
as examples. We also discuss, in Section HI the relationship of such results to constraints 
arising from other quantum-mechanical expressions, such as the famous Bethe sum rule 
[IS] and other interconnections between the S p (n). Finally, in Sec. [5j we briefly discuss 
mult i- summation expressions, generalizing Eqn. (jlj), motivated by constraints arising 
from higher orders of perturbation theory. 

2. The quantum bouncer and sum rule constraints 

The Schrodinger equation for the potential in Eqn. ([I]) is 

h 2 cPipJx) 

- 7: A3 V Fxt^Jx) = E n ip n (x) for < x < 00 (5) 

2m ax z 

where ip n {x) and E n are the energy eigenfunctions and eigenvalues, respectively. The 

appropriate boundary conditions are ip n (x = 0) = and ip n (x — ¥ 00) = 0. A change of 

variable, x = p(, transforms Eqn. (jSJ) into 

d ^ 1 = (C-UM0 (6) 
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where 



P 



(JL\ 



1/3 



and 



E n = 



'h 2 F 2 ' 



1/3 



C„ = £oC„ • (7) 



\2mF J " V 2m 

The solutions of Eqn. ([6]) are Airy functions, namely Ai(( — Cn) and -Bi(C — Cn), and 
only the Ai(Q solution is acceptable since the ip n (x) must be square integrable over 
the range (0, oo). The boundary condition at the infinite wall (C = 0) imposes the 
additional constraint that Ai(—( n ) = 0, so that the — ( n are simply the zeros of the 
relevant Airy function, — ( n . Since ( n = ( n , the energy eigenvalues are then directly 
given by E n = ( n £ and familiar WKB arguments or handbook results [13] can be used 
to derive an approximate formula for Cn for large n (quantum number), namely 



Cn 



Y(»-l/4) 



2/3 



(8) 



The normalization of these states, found first numerically in Ref . [3] , and then confirmed 
analytically in Refs. [5] and [8], make use of earlier published results (Refs. [6] and [7]), 
and can be written in the form 

Ai(( - Cn) 




MO 



(9) 



s/p\Ai'(-Q\ 

if we include the proper dimensional constant. 

Goodmanson |8j has found a recursion relationship for the position-matrix elements, 
which when written in terms of the normalized eigenstates and the scaled variable 
C = x/p, reads 



2<*i 



\n-k+l 



,P\ 



p{p - l)(p - 2)(p - 3)(n|C p - 4 |A;) + Ap(p - l)( ave (n\(?- 2 \k) 
- 2p(2p - lXnl?- 1 !*) + (Cn - tk) 2 (n\C P \k) 



(10) 



where Ca^e = (Cn + Cfc)/2. It is understood that for a given value of p, any expectation 
values of negative powers of C are to be ignored. 

Using this algorithm, one can find the first few diagonal and off-diagonal matrix 
elements needed for various sum rule calculations. For example, one has 



P 



P 



P 



P 



P 



2 : 



3: 



4: 



5: 



6: 



(n\(\n) = -y 



(n\( 2 \n 
(n\( 3 \n 

WC 4 |n 
(n\( 5 \n 



and 



P 



(n\(\k) 



8_C 
15 
16C 



35 

128C 
315 

256CJ 
693 



3 
7 



80Cn 

63 

1808C; 

3003 



\n— fc+1 



(Cfe _ Cn 



(11) 

(12) 
(13) 
(14) 
(15) 

(16) 
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94 (_1 \n-k+l 

p = 2: (n\( 2 \k) = { _ j (17) 



3|r,\ ( i \n— fc+1 



p = 3: (n|C d |A;) = (-1) 



720 48C„ 24 



p = 4: (n|C 4 |fc) = (-1 



,71— fe + 1 



(C^ - Cn) 6 (C fe - Cn) 4 (Cfc - Cn) 3 

40340 3840C„ 1920 



;is) 



.(19) 



.(a-Cn) 8 (Ck-Cn) 6 (a -co 5 . 

(We can then insert the appropriate powers of p as needed for dimensional correctness 
in physical matrix elements.) 

We note that Eqn. (fTTj) is consistent with the quantum-mechanical virial theorem 
as it gives 

(n\V(x)\n) = Fp{n\(\n) = S (^f\ = 2 -E n (20) 

which is known to be appropriate for a symmetric power-law potential, Vk(x) = Vo\x/a\ k , 
with k — 1. Other physically useful matrix elements, such as those involving p 2 , can be 
obtained by writing p 2 = 2m(H — Fx). For example, this allows for the evaluation of 
diagonal matrix elements such as 

(n\f\n) = 2mS ( ^ J and (n\p 4 \n) = (2m£ ) 2 ( ^ J (21) 

or off- diagonal ones such as 

Af (_-\\n-k+l 

(n\p 2 \k) = -2mFp(n\C\k} = ,°! \. 2 . (22) 

(Cfe - Cn) 

Using these results, we can begin cataloging the constraints which result from the 
application of various well-known quantum-mechanical sum rules. For example, using 
the energy differences E k — E n = £o(Cfc — Cn) and the off-diagonal dipole matrix elements 
in Eqn. f TTB]) . we start with the most-cited sum rule, the one formulated by Thomas- 
Reiche-Kuhn as in Eqn. ()3]). We find that all of the dimensional parameters cancel, 
leaving the simplest constraint we encounter, namely 

*<»>=£(o4H- (23) 

This example and all of the other constraint equations for the S p (n) derived here, as 

well as the on-diagonal and off-diagonal matrix elements in Eqns. ( TTT1) - ( fl9|) . can be 

verified to essentially arbitrary accuracy by using the numerical ability of mathematical 

manipulation programs such as Mathematica® . 

The so-called 'monopole sum rule,' which has been used in applications to nuclear 

collective excitations [13], is given by 

2h 2 
J2(E k - £„)|(n|x 2 |A;>| 2 = — (n\x 2 \n) , (24) 

T m 

and is of a similar form to the TRK sum rule. Using the ( 2 matrix elements from 
Eqns. (fT2l and ffT7|) . we find the constraint equation 

^Sw^- (25) 
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We note that the monopole sum rule in Eqn. 
form derived by Wang [T5j . 

2m 



is a special case of a more general 



XX£fc-£n)Kn|F(x)|A;)| : 

k 



n 



dF(x) dF* 



x 



dx 



dx 



n 



(26) 



which simplifies if the function is real so that F(x) = F*(x). This general result can 
also be used to immediately reproduce the TRK sum rule by using F(x) = x. 

Bethe and Jackiw [16j [17] derive several other sum rules for dipole moment matrix 
elements by using multiple commutation relations with the Hamiltonian, thus yielding 
higher powers of the energy difference. These higher-order sum rules include: 



z2(E k — E n 



{n\x 



h 2 

-^(n\p 2 \n) 



77T 



in 2 

— [E n 
m 



J2(E k - E n f\{n\x 



2m 2 



n 



d 2 V{x) 



and 



h 4 



Y,(E k - E n f\{n\x\k)\ 2 = — 



n 



dx 2 



'dV(x) 
dx 



n 



(n\V(x)\n)} (27) 



(28) 



?? 



(29) 



where Eqns. ( 128]) and (129]) are sometimes called the "force times momentum" and " force 
squared' sum rules, respectively We recall that not all such sum rules are guaranteed 
to lead to convergent expressions. 

The first of these three higher-order sum rules, Eqn. (1271) . gives the relation 

1 __Cn 

3 ' 



S-2 



k^n vVc 



>>n/ 



(30) 



This sum rule leads to a convergent result since, for large k, the terms in the summation 
scale as l/fc 4//3 . The sum rule results in Eqns. ( 128]) and ( 129]) . however, do not converge 
due to the derivatives of the discontinuous potential energy function in Eqn. ([[]). This 
implies that Si(n) is not convergent. 

As mentioned above, it has been emphasized that the standard expression for the 
second-order energy shift in perturbation theory due to an added potential energy term 
of the form V(x), given by 

\(n\V(x)\k}\ 2 



EP 



E 

k^n 



E& - E 



(0) 



(31) 



is also a form of energy-weighted sum rule. The authors of Ref. [12] have used this 
fact to evaluate the Stark shift, that is the second order energy shift due to an external 
constant field, with a potential of the form V(x) = Fx, in two model systems, the 
infinite square well and single attractive ^-function potential. Then, using the same 
mathematical techniques as for the confirmation of many other sum rules in those two 
cases (where the relevant tools are the Mittag-Lefler theorem and standard contour 
integration methods respectively) one can evaluate E^ in closed form. 
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In this situation, the addition of a uniform external field V(x) = Fx to the potential 
in Eqn. ([1]) leads to a soluble problem with a simple redefinition of the constant force, 
F — > F + F, giving an exact value for the new energy eigenvalues, 



E = C 

J - J n Sn 



[ h 2 (F + F) 2 



1m 



1/3 



E„ 



n 



^ 



2/3 



where E n = £ ( n . (32) 



The term in brackets can be easily expanded giving predictions for the first-, second-, 
and third-order perturbation theory results, namely 

E® = I (f ) (%) , ^ = ~(f) 2 (£ba, ^ 3) = ^(^) 3 ^Cn). (33) 

The first-order result is easily confirmed by noting that 

E& = (n\Fx\n) = Fp(n\(\n) = (|) (Fp) (^A = \ (|) (S Q (34) 




where we use the diagonal dipole matrix element in Eqn. (ITT]) . The second-order shift 
equation can then be used as a new constraint on a different combination of Airy function 
zeros, giving 

SJn) = Y - 1 -— = ^ (35) 

and we note that S% (n) = 11S^,(n). Higher-order perturbative corrections can in 
principle be used to derive closed form expressions for more complex combinations of 
inverse powers of ((k — Cn), since one can expand the exact result in Eqn. (1321) to 
arbitrarily high order. In fact, we use the third-order expression to briefly discuss 
multi-index summation generalizations of Eqn. (j4]) in Sec. EJ 

To briefly summarize, the evaluation of several well-known quantum-mechanical 
sum rules and the evaluation of the related energy-weighted sum over dipole matrix 
elements from perturbation theory have provided closed form expressions for S p (n) for 
p = 1,3,5 and 7. In the next section, we show how to systematically evaluate S p (n) for 
all natural p > 1. 

3. Systematic method for construction of S p (n) 

Before proceeding, we recall the techniques that are used in the derivation of many 
of the familiar quantum-mechanical sum rules, especially those involving dipole matrix 
elements. As an example, for the TRK sum rule in Eqn. ([3]), we write the standard x,p 
commutation relation, [p, x] = —in, bracketed by energy eigenstates: 

— ih = (n\px — xp\n) = J^ {(n\p\k) (k\x\n) — (n\x\k) (k\p\n)} (36) 

all k 

where we have also inserted a complete set of states to obtain the right-hand side of the 
equality. We can then use a second commutation relation, namely 

[H,x] = ^-[p 2 ,x] = —p, (37) 

1m mi 
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where we assume a standard one- dimensional Hamiltonian of the form H = p 2 /2m + 
V(x), to write 

/ i -in im i uii 11 n im(E n -E k ) 

(n\p\k) = — {n\[H, x\\k) = {n\x\k) . (38) 

There is a similar expression for (k\p\n) and combining these two results in Eqn. ( 1361) 
gives the TRK sum rule. Note that Eqn. (13"8"|) gives (n|p|A;) = if n = k as is appropriate 
for energy eigenstates where the average momentum should vanish in a stationary state. 
Specializing now to the case of the quantum bouncer, and using the result in 
Eqn. (fl6|) . we find the more specific result, 

(n\p\k) = - \/ ... (39) 

and we note that (n\p\k)* = (k\p\n). This is the lowest-order (in inverse powers of 
Cfc — Cn) term possible, and is our starting point. 

We see that the important ingredients are insertion of a complete set of states, 
the matrix-element connection in Eqn. (I38p . and appropriate commutation relations. 
Motivated by these methods, we start with the simple closure relationship for the 
momentum operator, namely 

^2(n\p\k) (k\p\n) = (n\p 2 \n) . (40) 

all k 

The only non-zero elements on the left-hand side are the off-diagonal matrix elements 
in Eqn. ( |39l) . while the diagonal matrix element on the right-hand side can be evaluated 
using Eqn (|2TJ) . Inserting these results gives 

S 2 (n) = | . (41) 

This is the same result as obtained from the sum rule in Eqn. (12TJ) . which is correct 
since that expression is most simply obtained from the closure relationship in Eqn. ( 140]) . 
using the expression in Eqn. (13"8l twice. 

The next higher power of (( k — Cn) -1 is obtained by inserting a complete of states 
into the commutation relation [x, p] = ih, namely 

y^ {(n\x\k) (k\p\n) — (n\p\k) (k\x\n)} = (n\ih\n) = ih (42) 

all k 

which, of course, reproduces the TRK sum rule result, giving S 3 (n) = 1/4. We note 
here that in any summation involving (n\p\k), the n = k term is not present. 
To evaluate S^n), we use the x-closure relationship, 

'^2(n\x\k)(k\x\n) = (n\x 2 \n) (43) 

all k 

and by explicitly including both on- and off-diagonal terms on the left-hand side (and 
removing all dimensional constants) this gives us 

MCH| 2 + E ^ \^ = (n\C 2 \n) (44) 
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or 



'f ) + 4S <<"> = § 



yielding 



c 

v ; 45 



(45) 



The correct iterative procedure required to evaluate S p (n) for any natural value 
of p is now clear. The off-diagonal matrix elements for (n\( q \k) will have a leading 
term of order ((j, — ( n )~ 2q . We then apply closure to the general commutator result 
[x q ,p] = iqhx q ~ l in the form 

y^ {(n\x q \k) (k\p\n) — (n\p\k) (k\x q \n)} = iqfi(n\x q ~ l \n) . (46) 

all k 

The recursive relation in Eqn. (ffOl) can then be used straightforwardly (if tediously) to 
evaluate both (n\x q \k) and (n\x q ~ l \n) to obtain a closed- form expression for 5*25+1 (n) m 
terms of explicit powers of ( n and values of S p (n) for p < 2q. This procedures allows one 
to increment the value of p by one, since the inclusion of the (n\p\k) matrix elements 
adds one more inverse power of ((k — ( n ). 

To evaluate S p (n) for p values two units higher, we generalize the x-closure 
relationship in Eqn. fT4*3|) to 

J2(n\x q \k)(k\x\n) = (n\x q+1 \n) (47) 

all k 

and again evaluation of only x-dependent on- and off-diagonal matrix elements will 
suffice to obtain an expression for S2 q +2{ n ), as the (k\x\n) term adds two inverse powers 
of (Cfe-Cn)- 

In this way, we are able to systematically evaluate S p (n) for all p > 1, since the 
p = 1 case, corresponding to use of Eqn. (128]) . does not converge. We present below 
results obtained in this way for the cases p 



2, ...,11, 



S 2 (n 

S 3 (n 

S±{n 

S 5 (n 

S 6 (n 

S 7 (n 

S 8 {n 

S 9 (n 

Sio(n) 

Su(n) 



+ 



3 
1 

4 

Vrt 

45 

(n 

36 

945 

c 2 

Sn 

270 

c 4 

4725 

C 3 
2100 

93555 

c 4 



112 



+ 



+ 



5Cn 

2268 
1 



+ 



+ 



2240 
6HC 2 



1496880 
43Cn 



17010 272160 



(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
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all of which have also been verified numerically. 

4. Other sum rule constraints 

There are a large number of other possible constraints, obtainable from both the use of 
the systematic approach developed above, or through other familiar sum rules. 

For increasing values of p, there are often several closure- motivated constraints that 
can be applied to obtain closed- form expressions for S p (n). For example, to evaluate 
Ss(n), we can use either 

'^2(n\x s \k)(k\x\n) = (n\x 4 \n) = ^2,(n\x 2 \k)(k\x 2 \n) (58) 

all k all k 

and use one as a cross-check against the other. Alternatively, such identities can be 
used as further constraints among the various S p (n). The expression above, for example, 
requires that 

{16£ , 3\ /2C„ X 



35 7 3 



+ 1440S 8 (n) - 96( n S e (n) - 47S 5 (n) 



= {jpj +576S 8 (n) (59) 

which is easily verified. 

One of the more famous constraints on dipole matrix elements is the Bethe sum 
rule [IE] 

J2(E k -E n )\(n\e^\k)\ 2 = ^, (60) 

k 

which was developed in the study of energy loss mechanisms, eventually leading to 
the Bethe-Bloch formula. This single sum rule actually provides an infinite tower of 
constraints on the S p (n). To see this, we expand the matrix element of the exponentials 
via 

2 3 4 

— (n\x 2 \k) — i—(n\x 3 \k) -\ — - 
2! 3! 4! 

The odd-order (in q) terms on the left-hand side of Eqn. (1601) automatically cancel since 

it is an even function by construction, while the O(q ) terms is absent since (n\k) = if 

n j^ k. The lowest non-vanishing term, the 0(q 2 ) term on the left-hand side, saturates 

the right-hand side by reproducing the TRK sum rule. All higher-order terms must 

therefore vanish. For example, the vanishing of the 0{q A ) term implies that 

\j2( E k~ E n )(n\x 2 \k)(k\x 2 \n) = i]T(£ fc - E n ) (n\x 3 \k) (k\x\n) (62) 

4 k d k 

which is a different constraint equation involving Sg(n) and including lower-order terms 
than the one in Eqn. (158]) . 



(n\e iqx \k) = (n\k) + iq(n\x\k) - \(n\x 2 \k) - i?-(n\x 3 \k) + \(n\x 4 \k) + ■■■ .(61) 
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Closure relationships involving mixed combinations of x and ft of higher order in the 
momentum operator are also possible, often reproducing earlier results. For example, 
the closure-relation for p A , namely 

'^2(n\p 2 \k)(k\p 2 \n) = (n\p 4 \n) (63) 

all k 

can be used in conjunction with Eqns. (I2T]) and (|22|) to evaluate 64(71). Wang [15] has 
also derived a number of sum rules involving mixed position- and momentum-matrix 
elements which can be used. 

5. Multi-index summations 

To explore more complex relationships involving inverse powers of (£& — ( n ), we extend 
the perturbation theory analysis of the Stark effect for the quantum bouncer to third 
order. We first recall that the third-order correction in perturbation theory due to a 
general V(x) term is given by 

(n\V(x)\k)(k\V{x)\j)(j\V(x)\n) 



4 3) =T,T, 



(°) Z7(°)VZ7(°) p(°)> 






\(n\V(x)\k)\ 



2 



(n|V(x)|n)E "(0) ^c ' ■ (64) 

fc^„ (£/n - % J 



Using V^ = Fx as above, the exact third-order result in Eqn. ( 1331 . and removing all 
dimensional factors gives the constraint 

fC = j, (n|C|fc)(fc|Cb>(7|C|n> | y. \(n\(\k)\ 2 (k\(\k) 

°* k^j^n (ik — (,n){{j — in) j=k^n [Qk ~ Qn) 

- (n\C\n) 2_s (r _ t \i - ( 65 ) 

k^n \^ k W 

In evaluating the double sum, we have separated off the cases where j = k ^ n, leaving 
the distinct j 7^ k 7^ n terms. 

Inserting the appropriate diagonal and off-diagonal matrix elements in Eqns. (TTTT) 
and (fl6"|) . we then find the relationship 

- 8 £. ft-w-w-u- + i |Ss(n) + C " S6( " )1 " ^ S6(,I) = « (66) 

or 

rsA3(n) = J^ (Cfc _ Cn)3(Cfc _ 0)2(c ._ Cn)3 = 3^ > (67) 

a remarkably simple identity, which we have also confirmed numerically. 

The presence of double (or higher) summations can be systematized, as in Sec. El 
by introducing more than one insertion of a complete set of states into any quantum 
identity, either resulting from simple closure or commutation relations. For example, 
we can write 

52 ^2( n \ x \k) (k\ x \j) (j\ x \ n ) = (^l^ 3 !^) • (68) 

all k all j 
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In evaluating the double sum, we must consider the following special cases, namely (i) 
n = k = j, (ii) two equal contributions arising from n = k ^ j and n = j 7^ k, (iii) 
j = k j^ n, and (iv) the completely distinct double sum where j 7^ k 7^ n. Including all 
of these possibilities gives the constraint 

+ 2 (^W) + I [S 3 (n) + CnS^n)} 




hn (C* " Cn) 2 (0 " G) 2 (Ci - Cn) 2 35 + 7 l > 



or 



2t 3 5 
T 2 2 2 = -2L + — (70) 

' ' 945 168 V ; 

which we have also confirmed numerically. Clearly an infinite number of multi-index 

summations can be evaluated in closed form in this way. 

6. Conclusions and future directions 

In conclusion, motivated by identities derived using a variety of quantum-mechanical 
sum rules, we have developed techniques to systematically evaluate sums of the form 
in Eqn. (J3J) for arbitrary natural p > 1, explicitly exhibiting results for S p (n) for 
p = 2, ...,11. In addition, we have identified many other additional constraints on 
the S p (n) arising from self-consistency of the method and the Bethe sum rule. Using 
higher-order perturbation theory results as a starting point, we have also discussed the 
existence of mult i- summation constraints arising from the the repeated use of insertion 
of a complete set of states. Future work might involve attempts at inductively generated 
closed form expressions for the S p (n) as well as exploration of the algebraic structures 
suggested by the pattern of results in Eqns. (l4"8"j) - (1571) . where each term only contains 
powers of ( n modulo three. 
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